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ABSTRACT 

We discuss implementation of TOPbase and Iron Project opacities for stellar spec- 
tral codes. We use a technique employed by Peach, where a Boltzmann-averaged cross 
section is calculated for selected temperatures, and the opacity obtained from double 
interpolation in temperature and wavelength. It is straightforward to include all lev- 
els for which cross sections have been calculated. Boltzmann-averaged cross sections 
for Fe II show a local maximum between 1700 and 2000A. We suggest this feature 
arises from id^Asnl — > Sd^Apnl transitions within Fe ii. lUE spectra of iron-rich CP 
stars show local minima in this region. Theoretical calculations of a representative 
stellar continuum demonstrate that Fe ii photoionization contributes significantly to 
the observed minima. 
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1 INTRODUCTION AND RATIONALE 

State of the art Photoionization cross sections have been 
available from The Opacity Project (Seaton 1994, The Opac- 
ity Project Team, 1995) for more than a decade. Neverthe- 
less, stellar atmospheres codes have been slow to fully in- 
corporate these results. Surely, one reason for the delay in 
adopting these calculations is that metal opacities are not 
particularly important for many stellar atmospheres. The 
well-known opacities from hydrogen, bound-free, and H~, 
usually dominate other sources. 

But it is also tedious to incorporate the detailed TOP- 
base results into workable subroutines with typical algo- 
rithms, which treat levels individually, or in groups. This 
technique is difficult for complex atoms, where no simple 
parameterization can describe all of the cross sections ade- 
quately. In this paper, we propose adoption of the method 
used by Peach (1967), which provides a practical solution to 
this problem. We then discuss an example of stellar contin- 
uum calculations that uses TOPbase/Iron Project cross sec- 
tions along with her algorithm. We first outline the general 
problem of dealing with the TOPbase/Iron Project material. 

Photoionization cross sections in TOPbase are given for 
typically hundreds of levels. For each level, the cross sections 
may be given several hundred to a thousand photon energies. 
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The cross sections include large deviations from the smooth 
hydrogenic functions, that are possible in multi-electron sys- 
tems. The resulting features may be quite sharp, being es- 
sentially spectral lines, or they can be broader, when there 
is strong interaction between bound and continuum states. 
Not only do Fano profiles appear, but also broad, deep ab- 
sorption minima (Cooper minima), which result from nodes 
in the radial wave functions. These are especially appar- 
ent in transitions involving s-electrons, as in Na I. Broad 
absorption maxima can result from photoexcitation of core 
electrons (PEC), essentially bound-bound transitions of in- 
ner electrons. These maxima can appear in the cross sections 
of numerous levels. 

A consensus has emerged that it is permissable to 
smooth over the resonances of photoionization cross sec- 
tions (cf. Bautista, Romano, and Pradhan 1998, Prieto et 
al. 2002). This smoothing affects the sharp features which 
we shall call high-frequency components, but leaves broader 
ones, such as Cooper minima or PEC. Prieto has published 
a number of resonance-smoothed photoionization cross sec- 
tions (RSP) on his website: 

|http:/ /hebe.as.utexa s.edu/at/at.cgi Our current opacity 
routines use the Prieto RSP's for all atoms and ions with 
the exception of Fe I and ll. Specifically, the following atoms 
and atomic ions are included: Li I, Be I and ll, B I and ll, C 
I, Na I, Mg I and ll, Al I, Si I and ll, S I, and Ca I and ll. 

For Fe I, we have sampled each lOOA from material dis- 
cussed by Bautista (1997). The Fe ll cross sections are from 
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a web site of Dr. Nahar 

|http://www-astronomy.mps.ohio-state.edu~nahar/ cf. Na- 
har and Pradhan 1994). 

While the RSP's omit some of the higher TOPbase lev- 
els, they may still contain several hundred levels (e.g. Si l), 
with several hundred data points for each level. If all of this 
information were included in a standard FORTRAN data 
statement, the result would be thousands of records long. 
However, for practical calculations, the individual cross sec- 
tions are not necessary. We need only their appropriately 
weighted sums. 



! 



Fell 

T = 7000K 
Dotted: Hydrogenic 



2 ALGORITHM AND CALCULATIONS 

The continuous opacity due to a given atom or atomic ion 
is the weighted sum of the cross sections over a number of 
levels. We define the K-factor: 



Kfact = ^ ffn ■ exp [-{xn/kT)] ■ (t„(A) 



(1) 



Here, g„ is the statistical weight of the level n, Xn/kT the 
appropriate Boltzmann factor, and a„{X) is the photoion- 
ization cross section. 

Consider the opacity due to photoionizations from level 
n of an atom or atomic ion in the i*'' ionization stage. Write 
the corresponding number density as Nn,i- Put Nt,i for the 
sum of the N„,i's in the i**^ stage of ionization. We have 



iVn.i = Nt.i 



9n 



■exp(-Xn/fcT), 



(2) 



where Ui{T) is the partition function for the ionization 
stage. The opacity due to all levels from this ion is thus 



■ Kfact. 



(3) 



In LTE, this K-factor is a function only of the temper- 
ature and wavelength. Since it is a sum over cross sections, 
some of the high-frequency structure for individual levels 
will be smoothed. One sees primarily, the absorption edges, 
and the most pronounced structure, for example, the PEC. 
The algorithm consists of computing this K-factor over the 
wavelengths likely to be or relevance, say from 500 to 8000A, 
for a relatively small number of temperatures. It is straight- 
forward to include all of the levels in TOPbase, or the RSP's. 

The opacities for any given atom or atomic ion are 
then found from bilinear interpolations in this temperature- 
wavelength grid. 

We have calculated the K-factor at 20A intervals, for 
all species except for Fe I. 

With 20A sampling, K-factors were first obtained for 
376 points [(8000-500) /20 -I- 1]. Then, unnecessary points 
were eliminated wherever "interior" data could be fit by lin- 
ear interpolations from bounding points to within 2%. For 
example, if the data at A4020 could be fit by interpolation 
between the points at AA4000 and 4040, the point at A4020 
was dropped. This technique resulted in considerable short- 
ening of the data statements, especially for lower temper- 
tures. 

The current routines for Fe I and ll are considered pro- 
visional, because they have not been resonance smoothed. 
The Nahar Fe ll photoionization file is some 33 megabytes 
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Figure 1. K-factors for Fe ll 



in size, and contains cross sections for 745 states (or "LS- 
multiplets"). We computed K-factor's every 20 A (AA500- 
8000) for 8 temperatures (2000 - 16000K). The adopted cross 
section at each of the grid wavelengths is an average of the 
nearest four TOPbase points. We have postponed improv- 
ing this algorithm until smoothed cross sections are avail- 
able. The current Fe ll subroutine is some 360 FORTRAN 
statements, including data statements, comments, and com- 
mented diagnostics. Comparisons between the K-factors for 
smoothed and unsmoothed cross sections for lighter species 
show differences in detail, but not major trends, extending 
over tens of angstroms. The feature discussed below is rela- 
tively broad. 

Fig.Qshows plots of the K-factors for Fe ll at 7000 and 
12000K, based on the Nahar data file. The dotted line shows 
the resulting "H-factors," that is, where the photoionization 
cross sections are all assumed hydrogenic (Cowan 1981, Eq. 
18.47). The feature between about 1700 and 2000A is no- 
table. We suggest that it arises from core photoexcitations 
3d^4snZ SdPipnl within Fe ll. It is important to notice 
that in these sort of transitions the outer nl electron acts as 
spectator and has only little effect on the wavelengths and 
radiative rates of the core transitions. The related 3d^4s 
3d^4p transitions in Fe ill give rise to the strong lines in 
UV Multiplets 34, 50, 51, and 52. We have made a further 
investigation of the nature of the Sd^Asnl Sd^ipnl tran- 
sitions in Fe ll using the atomic structure code AUTOSTRUC- 
TURE (Badnell 1986, 1997; Badnell & Pindzola 1989). This 
code is an extension of the atomic structure program super- 
structure by (Eissner et al. 1974). It allows for the calcu- 
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CALCULATED FLUXES 
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Figure 3. Intensity factors for Fell lines subject to autoioniza- 
tlon. The ordinate, log(g/) — 0.5 ■ Xi is roughly proportional to 
the logarithm of the strength of an absorption line if saturation 
did not occur. The logarithmic Boltzmann factor, corresponds to 
e = 5040/T = 0.5 or 1080K. 
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Figure 5. Calculated continuum fluxes for Te = 11600K, and 
\og{g) = 3. The plots, from upper to lower show results for iron 
and silicon 0.1-, 1.0—, and 10.0-times solar abundances 



3 OBSERVATIONS 

lUE spectra of iron-rich stars show significant depressions in 
the wavelength range from roug hly 1700-1800A, as is shown 
in the top 4 panels of Fig. 2] These spectra were downloaded 
from the url http://archive.stsci.edu/iue/ Because of the 
high noise level of the spectra, we have chosen to plot only 
the maxima within lOA regions. Complete spectra may be 
seen at the url cited in the figure. Smith and Dworetsky 
(1993) report high iron abundances in 112 Her, and HR 
6000. The magnetic Ap star a^CVn has long been known to 
be iron rich (cf. Cohen 1970), while Muthsam and Cowley 
(1984) showed HR 6870 (HD 168733) to be quite iron rich. 
The stars 87 Psc and 53 Tau are iron poor (see Smith and 
Dworetsky 1993), and their spectra lack the strong minima 
seen in the other four stars. 

Lanz, et al. (1996) discuss the region on the short wave- 
length side of I6OOA. The interpretation is not straightfor- 
ward, though the figures shown here are clearly relevant for 
any interpretation. Note that HR 6000, notorious for its lack 
of strong Sill lines in ground-based photographic spectra, 
also lacks the depression short of A1600 (cf. Andersen, et al. 
1984). 



lation of level energies, and radiative and Auger rates within 
a statistical Thomas-Fermi-Dirac model potential (Eissner 
& Nussbaumer 1969). Thus, we compute the total of 8988 
L5'-coupling dipole allowed transitions Sd^isnl — > Scfiipnl 
for 5 ^ n ^ 10 and ^ I ^ 3, with the majority of the 
3d^4pnl states lying above the the first ionization thresh- 
old. A little less than half of these transitions, 4090, end up 
in upper states that can autoionize to Felll target states. 
Such transitions manifest themselves as PEC resonances in 
the close coupling photoionization cross sections of Nahar 
and Pradhan (1994). Several of the cross sections that con- 
tribute to the local peak in the 1700 to 2000A region are 
shown in Fig. |21 (following page). The remaining 4898 transi- 
tions go to levels that do not autoionize and must be treated 
as bound-bound transitions. We show a subset of the former 
lines in Fig|^ The concentration of transitions in the region 
I7OO-2OOOA is clear. 



4 APPLICATION TO FE II 

Fig. 13 shows calculated spectra in the wavelength region 
under consideration. The model was based on the adopted 
parameters of Kochukhov et al. (2002) , an effective temper- 
ature of 11600K and log(g) — 3.0. We also followed these 
authors in using crude abundances to set the model. We 
took abundances 0.1 SAD (Grevesse and Sauval 1998) for 
elements helium through neon, and 10 SAD for all heav- 
ier elements. This was only for the determination of the 
temperature-pressure structure. A flux-constant model was 
derived from an Atlas 9 model (Kurucz 1993) , but the depth 
depencence was recalculated using Michigan routines as de- 
scribed in earlier papers (cf. Cowley 1996), except that our 
opacity routines now incorporate TOPbase data for the 
species enumerated above. 

The three curves were all based on the same T{t), 
though a new pressure structure was calculated based on the 
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Figure 2. Fc II cross sections (Naliar and Pradlian 1996) for levels with wide maxima in the region I7OO-2OOOA. Hydrogenic cross 
sections are shown by the dotted curves. The first three digits at the top of each plot identify the level or term by multiplicity, L-value, 
parity (O=even,l=odd). The fourth digit is a sequential number, from low to high excitation, of the term described by the first three 
digits. For example, 4 3 4 describes an even ''F term, that is the fourth such term in order of excitation. In TOPbase, the first three 
digits are called ISLP and the fourth, ILV. Chi(n) is he excitation energy of the level in eV. 



different iron and silicon abundances. Thus, the flux is not 
strictly constant in the models without the enhanced abun- 
dances. We do not consider this significant in the present ex- 
ploratory calculations. The three curves shown are for iron 
and silicon enhanced by a factor of 10 (lower), solar (mid- 
dle), and depleted by a factor of 10 (upper). All elements 
other than iron and silicon were assumed to have the same 
abundances used in the AtlasO model. 

The calculations show clearly that the Fell continuum 
is significant in this region. Of the three calculated minima, 
the one at the shortest wavelength (about A1640) is not well 
matched to the observations. It is certainly Fe 11 and not Si, 
as a separate calculation (not shown) has revealed. There is 
a possibly corresponding feature to be seen in the spectra 
of a'^CVn and HR 6000, though much weaker than in the 
calculations. This region deserves further study. 

The two longer-wavelength minima, at approximately 



AA1700 and 1740 fit the observations of the iron-rich spectra 
reasonably well. There is a small discrepancy in the wave- 
lengths; the calculated minima are some 20 to 30A short of 
the stellar minima. But these discrepancies are 2% or less, 
well within the accuracy with which energy levels in Fell 
were calculated (cf. Nahar and Pradhan, op. cit. Table 1). 



5 CONCLUSIONS 

Some of the absorption in the region AA1700-2000 of the 
spectra shown in Fig. 0| will result from the strong transi- 
tions in Fell as well as the multiplets mentioned above in 
Felll. But we show here that a measurable contribution to 
the depression must be due to the specific shapes of the 
Fell photoionization cross sections. Note that the lowest of 
the levels shown in Fig.|5| at 9.38eV is the 76th level above 
ground. Any attempt to treat higher levels in Fe 11 by a hy- 
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Figure 4. The top four panels are lUE spectra of iron-rich stars. The lower two panels show stars with iron deficiencies. Only the highest 
points within lOA intervals are plotted. The original spectra are available on line at the site: http:/ /archive. stsci.edu/iue/j 



drogenic approximation or simple parameterization would 
miss the relevant structure. 

A number of loose ends need to be tidied before our 
metal opacity routines can be considered satisfactory. We 
need smoothed cross sections for Fel and ll, as well as other 
heavy species, such as Mn I and ll or Co I and ll that can be 
relevant for special stars. It is also necessary to make sure 
that the partition functions used with the TOPbase/Iron 
Project cross sections are compatible with those levels. Pro- 
visional, spot checks for a few spectra show that this is not 



likely to be a serious problem, as the calculated levels are 
usually within a few per cent of the laboratory values. Line 
opacity was not relevant for the current calculations. Even- 
tually, we must consider the overlap between the TOPbase 
"resonances" and predicted lines in the Kurucz data base. 



These are matters for future contributions. 
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